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' Abstract We study the Fermi level structure of 2 + 1 -dimensional strongly interact- 

^^ i ing electron systems in external magnetic field using the AdS/CFT correspondence. 

The gravity dual of a finite density fermion system is a Dirac field in the background 
of the dyonic AdS-Reissner-Nordstrom black hole. In the probe limit the magnetic 
system can be reduced to the non-magnetic one, with Landau-quantized momenta 
and rescaled thermodynamical variables. We find that at strong enough magnetic 
fields, the Fermi surface vanishes and the quasiparticle is lost either through a 
pH , crossover to conformal regime or through a phase transition to an unstable Fermi 

surface. In the latter case, the vanishing Fermi velocity at the critical magnetic field 
triggers the non-Fermi liquid regime with unstable quasiparticles and a change in 
^^ ' transport properties of the system. We associate it with a metal-"strange metal" 

1/^ , phase transition. We compute the DC Hall and longitudinal conductivities using the 

CO ' gravity-dressed fermion propagators. As expected, the Hall conductivity is quan- 

OO , tized according to integer Quantum Hall Effect (QHE) at weak magnetic fields. 

At strong magnetic fields, new plateaus typical for the fractional QHE appear Our 
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■"sj" , pattern closely resembles the experimental results on graphite which are described 
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using the fractional filling factor proposed by Halperin. 
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1 Introduction 

The study of strongly interacting fermionic systems at finite density and temperature 
is a challenging task in condensed matter and high energy physics. Analytical meth- 
ods are limited or not available for strongly coupled systems, and numerical sim- 
ulation of fermions at finite density breaks down because of the sign problem (l]. 
There has been an increased activity in describing finite density fermionic matter by 
a gravity dual using the holographic AdS/CFT correspondence |[2|. The gravitational 
solution dual to the finite chemical potential system is the electrically charged AdS- 
Reissner-Nordstrom (RN) black hole, which provides a background where only the 
metric and Maxwell fields are nontrivial and all matter fields vanish. In the classi- 
cal gravity limit, the decoupling of the Einstein-Maxwell sector holds and leads to 
universal results, which is an appealing feature of applied holography. Indeed, the 
celebrated result for the ratio of the shear viscosity over the entropy density lO is 
identical for many strongly interacting theories and has been considered a robust 
prediction of the AdS/CFT correspondence. 

However, an extremal black hole alone is not enough to describe finite density 
systems as it does not source the matter fields. In holography, at leading order, the 
Fermi surfaces are not evident in the gravitational geometry, but can only be detected 
by external probes; either probe D-branes ||2] or probe bulk fermions lH |5] IS E) . 
Here we shall consider the latter option, where the free Dirac field in the bulk carries 
a finite charge density 18J. We ignore electromagnetic and gravitational backreac- 
tion of the charged fermions on the bulk spacetime geometry (probe approximation). 
At large temperatures, T ^ ji, this approach provides a reliable hydrodynamic de- 
scription of transport at a quantum criticality (in the vicinity of superfluid-insulator 
transition) ID. At small temperatures, T -^ jj., in some cases sharp Fermi surfaces 
emerge with either conventional Fermi-liquid scaling fF] or of a non-Fermi liquid 
type ID with scaling properties that differ significantly from those predicted by the 
Landau Fermi liquid theory. The non-trivial scaling behavior of these non-Fermi 
liquids has been studied semi-analytically in fJ] and is of great interest as high-r^ 
superconductors and metals near the critical point are believed to represent non- 
Fermi liquids. 

What we shall study is the effects of magnetic field on the holographic fermions. 
A magnetic field is a probe of finite density matter at low temperatures, where the 
Landau level physics reveals the Fermi level structure. The gravity dual system is 
described by a AdS dyonic black hole with electric and magnetic charges Q and 
H, respectively, corresponding to a 2 + 1 -dimensional field theory at finite chemical 
potential in an external magnetic field ifTOl . Probe fermions in the background of the 
dyonic black hole have been considered in ifTTl [T2l ; and probe bosons in the same 
background have been studied in il3J. Quantum magnetism is considered in llT4l . 

The Landau quantization of momenta due to the magnetic field found there, 
shows again that the AdS/CFT correspondence has a powerful capacity to unveil 
that certain quantum properties known from quantum gases have a much more ubiq- 
uitous status than could be anticipated theoretically. A first highlight is the demon- 
stration ifTSl that the Fermi surface of the Fermi gas extends way beyond the realms 
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of its perturbative extension in the form of the Fermi-liquid. In AdS/CFT it appears 
to be gravitationally encoded in the matching along the scaling direction between 
the 'bare' Dirac waves falling in from the 'UV boundary, and the true IR excitations 
living near the black hole horizon. This IR physics can insist on the disappearance 
of the quasiparticle but, if so, this ' critical Fermi-liquid' is still organized ' around' a 
Fermi surface. The Landau quantization, the organization of quantum gaseous mat- 
ter in quantized energy bands (Landau levels) in a system of two space dimensions 
pierced by a magnetic field oriented in the orthogonal spatial direction, is a sec- 
ond such quantum gas property. We shall describe here following ifTTI . that despite 
the strong interactions in the system, the holographic computation reveals the same 
strict Landau-level quantization. Arguably, it is the mean-field nature imposed by 
large A^ limit inherent in AdS/CFT that explains this. The system is effectively non- 
interacting to first order in 1 /N. The Landau quantization is not manifest from the 
geometry, but as we show this statement is straightforwardly encoded in the sym- 
metry correspondences associated with the conformal compactification of AdS on 
its flat boundary (i. e., in the UV CFT). 

An interesting novel feature in strongly coupled systems arises from the fact that 
the background geometry is only sensitive to the total energy density Q^ +H^ con- 
tained in the electric and magnetic fields sourced by the dyonic black hole. Dialing 
up the magnetic field is effectively similar to a process where the dyonic black hole 
loses its electric charge. At the same time, the fermionic probe with charge q is 
essentially only sensitive to the Coulomb interaction gqQ. As shown in ifTTI . one 
can therefore map a magnetic to a non-magnetic system with rescaled parameters 
(chemical potential, fermion charge) and same symmetries and equations of motion, 
as long as the Reissner-Nordstrom geometry is kept. 

Translated to more experiment-compatible language, the above magnetic-electric 
mapping means that the spectral functions at nonzero magnetic field h are identical 
to the spectral function at /i ~ for a reduced value of the coupling constant (fermion 
charge) q, provided the probe fermion is in a Landau level eigenstate. A striking 
consequence is that the spectrum shows conformal invariance for arbitrarily high 
magnetic fields, as long as the system is at negligible to zero density. Specifically, 
a detailed analysis of the fermion spectral functions reveals that at strong magnetic 
fields the Fermi level structure changes qualitatively. There exists a critical magnetic 
field at which the Fermi velocity vanishes. Ignoring the Landau level quantization, 
we show that this corresponds to an effective tuning of the system from a regular 
Fermi liquid phase with linear dispersion and stable quasiparticles to a non-Fermi 
liquid liquid with fractional power law dispersion and unstable excitations. This 
phenomenon can be interpreted as a transition from metallic phase to a "strange 
metal" at the critical magnetic field and corresponds to the change of the infrared 
conformal dimension from v > 1/2 to v < 1/2 while the Fermi momentum stays 
nonzero and the Fermi surface survives. Increasing the magnetic field further, this 
transition is followed by a "strange-metal"-conformal crossover and eventually, for 
very strong fields, the system always has near-conformal behavior where kf =0 and 
the Fermi surface disappears. 
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For some Fermi surfaces, this surprising metal-"strange metal" transition is not 
physically relevant as the system prefers to directly enter the conformal phase. 
Whether a fine tuned system exists that does show a quantum critical phase transi- 
tion from a FL to a non-FL is determined by a Diophantine equation for the Landau 
quantized Fermi momentum as a function of the magnetic field. Perhaps these are 
connected to the magnetically driven phase transition found in AdS5/CFT4 lfT6l . We 
leave this subject for further work. 

Overall, the findings of Landau quantization and "discharge" of the Fermi surface 
are in line with the expectations: both phenomena have been found in a vast array of 
systems iflTl and are almost tautologically tied to the notion of a Fermi surface in a 
magnetic field. Thus we regard them also as a sanity check of the whole bottom-up 
approach of fermionic AdS/CFT ||4]|6]|5][l5], giving further credit to the holographic 
Fermi surfaces as having to do with the real world. 

Next we use the information of magnetic effects the Fermi surfaces extracted 
from holography to calculate the quantum Hall and longitudinal conductivities. Gen- 
erally speaking, it is difficult to calculate conductivity holographically beyond the 
Einstein-Maxwell sector, and extract the contribution of holographic fermions. In 
the semiclassical approximation, one-loop corrections in the bulk setup involving 
charged fermions have been calculated |15|. In another approach, the backreaction 
of charged fermions on the gravity-Maxwell sector has been taken into account and 
incorporated in calculations of the electric conductivity IB]. We calculate the one- 
loop contribution on the CFT side, which is equivalent to the holographic one-loop 
calculations as long as vertex corrections do not modify physical dependencies of 
interest lITSlfTSl . As we dial the magnetic field, the Hall plateau transition happens 
when the Fermi surface moves through a Landau level. One can think of a differ- 
ence between the Fermi energy and the energy of the Landau level as a gap, which 
vanishes at the transition point and the 2 + 1 -dimensional theory becomes scale in- 
variant. In the holographic D3-D7 brane model of the quantum Hall effect, plateau 
transition occurs as D-branes move through one another [EOl . In the same model, a 
dissipation process has been observed as D-branes fall through the horizon of the 
black hole geometry, that is associated with the quantum Hall insulator transition. 
In the holographic fermion liquid setting, dissipation is present through interaction 
of fermions with the horizon of the black hole. We have also used the analysis of the 
conductivities to learn more about the metal-strange metal phase transition as well 
as the crossover back to the conformal regime at high magnetic fields. 

We conclude with the remark that the findings summarized above are in fact 
somewhat puzzling when contrasted to the conventional picture of quantum Hall 
physics. It is usually stated that the quantum Hall effect requires three key ingre- 
dients: Landau quantization, quenched disorderly and (spatial) boundaries, i. e., a 
finite-sized sample |2li. The first brings about the quantization of conductivity, the 
second prevents the states from spilling between the Landau levels ensuring the ex- 
istence of a gap and the last one in fact allows the charge transport to happen (as it is 



' Quenched disorder means that the dynamics of the impurities is "frozen", i. e. they can be re- 
garded as having infinite mass. When coupled to the Fernii liquid, they ensure that below some 
scale the system behaves as if consisting of non-interacting quasiparticles only. 
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the boundary states that actually conduct). In our model, only the first condition is 
satisfied. The second is put by hand by assuming that the gap is automatically pre- 
served, i. e. that there is no mixing between the Landau levels. There is, however, no 
physical explanation as to how the boundary states are implicitly taken into account 
by AdS/CFT. 

We outline the holographic setting of the dyonic black hole geometry and bulk 
fermions in the section |2] In section |3] we prove the conservation of conformal sym- 
metry in the presence of the magnetic fields. Section|4]is devoted to the holographic 
fermion liquid, where we obtain the Landau level quantization, followed by a de- 
tailed study of the Fermi surface properties at zero temperature in section |5] We 
calculate the DC conductivities in section|6l and compare the results with available 
data in graphene. 



2 Holographic fermions in a dyonic black hole 

We first describe the holographic setup with the dyonic black hole, and the dynamics 
of Dirac fermions in this background. In this paper, we exclusively work in the probe 
limit, i. e., in the limit of large fermion charge q. 



2.1 Dyonic black hole 

We consider the gravity dual of 3-dimensional conformal field theory (CFT) with 
global U{\) symmetry. At finite charge density and in the presence of magnetic field, 
the system can be described by a dyonic black hole in 4-dimensional anti-de Sitter 
space-time, AdS/i,, with the current J^ in the CFT mapped to a [/(I) gauge field Am 
mAdS. We use /i, v,p,... for the spacetime indices in the CFT andM,A^, . . . for the 
global spacetime indices in AdS. 

The action for a vector field Am coupled to AdSi, gravity can be written as 

5. = ^/^^xV^(^ + ^-|ft..^"^), (1) 

where g^p is an effective dimensionless gauge coupling and R is the curvature radius 
of AdS^. The equations of motion following from eq. ^ are solved by the geometry 
corresponding to a dyonic black hole, having both electric and magnetic charge: 

f2 d2 j 2 

ds^ = gMNdx^dx^ = -^ {-fdt^ + dx^ + dy^) + -j — . (2) 

The redshift factor / and the vector field Am reflect the fact that the system is at a 
finite charge density and in an external magnetic field: 



6 E.Gubankova, J. Brill, M. Cubrovic, K. Schalm, P. Schijven and J. Zaanen 

A, = M (l - -) , Ay = hx, A,- = A, = 0, (3) 

where Q and H are the electric and magnetic charge of the black hole, respectively. 
Here we chose the Landau gauge; the black hole chemical potential jx and the mag- 
netic field h are given by 

8fQ j^^EfH 

with ro is the horizon radius determined by the largest positive root of the redshift 
factor /(ro) = 0: 

M = rU^±^. (5) 

The boundary of the AdS is reached for r — > oo. The geometry described by eqs. 
(|2]|3]l describes the boundary theory at finite density, i. e., a system in a charged 
medium at the chemical potential fi ~ /Xbh and in transverse magnetic field h = h\±, 
with charge, energy, and entropy densities given, respectively, by 

The temperature of the system is identified with the Hawking temperature of the 
blackhole, rH-|/'(ro)|/4n:, 

Since Q and H have dimensions of [L]^, it is convenient to parametrize them as 

e' = 3r4, Q2+H^^3rl. (8) 

In terms of tq, r* and r„ the above expressions become 

f = ,+'-!»-i±^, (9) 



with 



r: 



2 



\/^* 



^ = V38F-^^, h^V38F ^ *^, * . (10) 

The expressions for the charge, energy and entropy densities, as well as for the 
temperature are simplified as 

_2V3_rl_ 1 rl+3rt/ro _ 27r r| 
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r = Az:ifi_4y (11) 

4;r R^\r^J ' ' 

In the zero temperature limit, i. e., for an extremal black hole, we have 

7^ = -^ ro = r„, (12) 

which in the original variables reads Q~ +H~ = 3rg. In the zero temperature limit 
(fT2] i. the redshift factor / as given by eq. (|9]l develops a double zero at the horizon; 

/^6 ^'~'"^ +^((r-r„)^). (13) 

As a result, near the horizon the AdS4 metric reduces to AdS2 x R^ with the curvature 
radius ofAdS2 given by 

«,.-L«^ ,14, 

This is a very important property of the metric, which considerably simplifies the 
calculations, in particular in the magnetic field. 

In order to scale away the AdS4 radius R and the horizon radius ro, we introduce 
dimensionless variables 

r -^ ror, r* -^ ror*, r„ -^ ror„, 

M^rlM, Q^rlQ, ,H^rlH, (15) 

and 

(f,x)^ — (f,x), Am^^Am, (o^^(o, 
ro R~ R 

2 

li^^li, h^§h, T^^T, 

ds-^R^ds^. (16) 

Note that the scaling factors in the above equation that describes the quantities of 
the boundary field theory involve the curvature radius of Aii5'4, not Arf5'2. 
In the new variables we have 



At = ^i(l--], jX = sj^gprl = gpQ, h = gpH, (17) 

and the metric is given by 



1 dr^ 



^'-■" -^■-^-^^^, ^ 



ds'-^r'-[-fdt^ + dx'- + dy'-) + ——^, (18) 
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with the horizon at r = 1, and the conformal boundary at r ^- oo. 

At r = 0, r** becomes unity, and the redshift factor develops the double zero 
near the horizon, 

(^:-i)V±2^ (19, 

n 
As mentioned before, due to this fact the metric near the horizon reduces to AdS2 x 
R^ where the analytical calculations are possible for small frequencies I?). However, 
in the chiral limit m = 0, analytical calculations are also possible in the bulk AdS4 
II22II . which we utilize in this paper. 



2.2 Holographic fermions 

To include the bulk fermions, we consider a spinor field \j/ in the AdS4 of charge q 
and mass m, which is dual to an operator ^ in the boundary CFT^ of charge q and 

dimension 

3 
A = -+mR, (20) 

with inR > —J and in dimensionless units corresponds to A = j+ m. In the black 
hole geometry, eq. (|2]l, the quadratic action for xj/ reads as 

S^ = i fd\x,/^{\irr'^&MW-mWW), (21) 

where i// = Xj/^'F-, and 

S!M = dM + - (OabMr"'' - iqAu, (22) 

where COatM is the spin connection, and F"'' = j[F",F^]. Here, M and a,b denote 
the bulk space-time and tangent space indices respectively, while /i,v are indices 
along the boundary directions, i. e. M = (r, /x). Gamma matrix basis (Minkowski 
signature) is given in fT\. 

We will be interested in spectra and response functions of the boundary fermions 
in the presence of magnetic field. This requires solving the Dirac equation in the 
bulk IE 13: 

{r'^&M-m)\ir^Q. (23) 

From the solution of the Dirac equation at small O), an analytic expression for the 
retarded fermion Green's function of the boundary CFT at zero magnetic field has 
been obtained in IJI. Near the Fermi surface it reads as Q: 

(0 — Vfki^ — L[(o,T) 
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where k± ^ k — kf is the perpendicular distance from the Fermi surface in mo- 
mentum space, hi and vp are real constants calculated below, and the self-energy 

Z = Ei+ il,2 is given by Q 

where v is the zero temperature conformal dimension at the Fermi momentum, 
V = Vkf, given by eq. ( fSSl ). pLq = jxq, /12 is a positive constant and the phase is 
such that the poles of the Green's function are located in the lower half of the com- 
plex frequency plane. These poles correspond to quasinormal modes of the Dirac 
Equation |23] and they can be found numerically solving F{(X)^,) = |i23|, with 

Fico)^ ^i — '^^^"-'-i^-^r (26) 

The solution gives the full motion of the quasinormal poles co* (kj^) in the complex 
(0 plane as a function of kj^. It has been found in 11231 171. that, if the charge of the 
fermion is large enough compared to its mass, the pole closest to the real (O axis 
bounces off the axis at ^x = (and co = 0). Such behavior is identified with the 
existence of the Fermi momentum kp indicative of an underlying strongly coupled 
Fermi surface. 

At r = 0, the self-energy becomes T^^g{(o/T) — > qO^*', and the Green's func- 
tion obtained from the solution to the Dirac equation reads |l7l 

where k±^ ~ \W- — kf. The last term is determined by the IR AdS2 physics near the 
horizon. Other terms are determined by the UV physics of the AdS^ bulk. 

The solutions to ( [23] l have been studied in detail in ||6l |5] [T] . Here we simply 
summarize the novel aspects due to the background magnetic field |fT9l 

• The background magnetic field h introduces a discretization of the momentum: 



k^k^s^ y/2\qh\l, with / £ A^, (28) 

with Landau level index I 12311121 . These discrete values of ^ are the analogue of 
the well-known Landau levels that occur in magnetic systems. 
There exists a (non-invertible) mapping on the level of Green's functions, from 
the magnetic system to the non-magnetic one by sending 



//2 



{H,Q,q)^{0,V^+H^,qJl--^^^j^_\. (29) 
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The Green's functions in a magnetic system are thus equivalent to those in the 
absence of magnetic fields. To better appreciate that, we reformulate eq. ( [29] l in 
terms of the boundary quantities: 

{h,^„T)^(^0,^„T(^l--^^y (30) 

where we used dimensionless variables defined in eqs. (I15I17| |. The magnetic 
field thus effectively decreases the coupling constant q and increases the chemi- 
cal potential jU = gfQ such that the combination jXij = jiq is preserved IITTI . This 
is an important point as the equations of motion actually only depend on this 
combination and not on jx and q separately ifTTl . In other words, eq. (|30t implies 
that the additional scale brought about by the magnetic field can be understood 
as changing jU and T independently in the effective non-magnetic system instead 
of only tuning the ratio jJ^/T. This point is important when considering the ther- 
modynamics. 

• The discrete momentum k^ff = y^2\qh\l must be held fixed in the transformation 
( |29] l. The bulk-boundary relation is particularly simple in this case, as the Landau 
levels can readily be seen in the bulk solution, only to remain identical in the 
boundary theory. 

• Similar to the non-magnetic system ifTTI . the IR physics is controlled by the 
near horizon AdS2 x R^ geometry, which indicates the existence of an IR CFT, 
characterized by operators ffi, I £ N with operator dimensions 5 = 1/2 + V/: 




v, = „,M„^ + lt^V^, (31) 



in dimensionless notation, and jx^ = jiq. AtT = 0, when r** = 1, it becomes 



V/ = i ^6{m^ + 2\qh\l)^^il (32) 

The Green's function for these operators ^/ is found to be ^/^(o) ^ aP'^' and 
the exponents V/ determines the dispersion properties of the quasiparticle excita- 
tions. For V > 1/2 the system has a stable quasiparticle and a linear dispersion, 
whereas for V < 1/2 one has a non-Fermi liquid with power-law dispersion and 
an unstable quasiparticle. 



3 Magnetic fields and conformal invariance 

Despite the fact that a magnetic field introduces a scale, in the absence of a chem- 
ical potential, all spectral functions are essentially still determined by conformal 
symmetry. To show this, we need to establish certain properties of the near-horizon 
geometry of a Reissner-Nordstrom black hole. This leads to the AdS2 perspective 
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that was developed in Q . The result relies on the conformal algebra and its relation 
to the magnetic group, from the viewpoint of the infrared CFT that was studied in 
ITTI . Later on we will see that the insensitivity to the magnetic field also carries over 
to AdS4 and the UV CFT in some respects. To simplify the derivations, we consider 
the case T = 0. 



3.1 The near-horizon limit andDirac equation in AdSi 

It was established in fT] that an electrically charged extremal Ac/^-Reissner-Nordstrom 
black hole has an AdS2 throat in the inner bulk region. This conclusion carries over 
to the magnetic case with some minor differences. We will now give a quick deriva- 
tion of the AdSo formalism for a dyonic black hole, referring the reader to |21 for 
more details (that remain largely unchanged in the magnetic field). 

Near the horizon r~r„ of the black hole described by the metric (2), the redshift 
factor f{r) develops a double zero: 

/W-6 ^'~["^ +e{{r-r„f). (33) 

Now consider the scaUng limit 

r~r„^X-^, r^A^'r, A ^ with T, ^finite. (34) 

In this Umit, the metric (|2]i and the gauge field reduce to 

ds^ = || {-dx^ + ^C') + ^ {dx^ + dy^) 

A,^ ^—f^y, A,^Hx (35) 

where R2 = -jr- The geometry described by this metric is indeed A£/52 x R'^. Phys- 
ically, the scaling limit given in eq. ( l34l i with finite T corresponds to the long time 
limit of the original time coordinate t, which translates to the low frequency limit of 

the boundary theory: 

a 

-^0, (36) 

where CO is the frequency conjugate to t. (One can think of A as being the frequency 
(O). Near the AdS^ horizon, we expect the AdS2 region of an extremal dyonic black 
hole to have a CFTi dual. We refer to JT) for an account of this Aii5'2/CFTi duality. 
The horizon of AdSi region is at ^ — > oo (coefficient in front of dz vanishes at the 
horizon in eq. ( |35] |) and the infrared CFT (IR CFT) lives at the AdS2 boundary at 
^ = 0. The scaling picture given by eqs. (1341135b suggests that in the low frequency 
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limit, the 2-dimensional boundary theory is described by this IR CFT (which is a 
CFTi). The Green's function for the operator G in the boundary theory is obtained 
through a small frequency expansion and a matching procedure between the two 
different regions (inner and outer) along the radial direction, and can be expressed 
through the Green's function of the IR CFT I?]. 

The explicit form for the Dirac equation in the magnetic field is of little interest 
for the analytical results that follow. It can be found in lfT9l . Of primary interest is 
its limit in the IR region with metric given by eq. (l35l i: 



(37) 

where the effective momentum of the Z-th Landau level is A/ = y^2\qh\l, fig = fiq 
and we omit the index of the spinor field. To obtain eq. (IJTI i. it is convenient to 
pick the gamma matrix basis as F? = — (J3, F* = /Ci and F' = —02. We can write 
expUcitly; 







(38) 



Note that the AdS2 radius R2 enters for the (t, Q directions. At the AdSx boundary, 
1^ — > 0, the Dirac equation to the leading order is given by 



Mi7^2''0 , R 

— 2 




The solution to this equation is given by the scaling function F'^'' ~ Ae+^ ^> 
Be_ ^^' where e± are the real eigenvectors of U and the exponent is 




v,^^J6{m^ + ^2\qh\l]R^-t^^. (40) 



The conformal dimension of the operator G in the IR CFT is 5i — j + Vi . Com- 
paring eq. ( l40l i to the expression for the scaling exponent in Q, we conclude that 
the scaling properties and the AdS2 construction are unmodified by the magnetic 
field, except that the scaling exponents are now fixed by the Landau quantization. 
This "quantization rule" was already exploited in ||231 to study de Haas-van Alphen 
oscillations. 
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4 Spectral functions 

In this section we will explore some of the properties of the spectral function, in 
both plane wave and Landau level basis. We first consider some characteristic cases 
in the plane wave basis and make connection with the ARPES measurements. 



4.1 Relating to the ARPES measurements 



In reality, ARPES measurements cannot be performed in magnetic fields so the 
holographic approach, allowing a direct insight into the propagator structure and the 
spectral function, is especially helpful. This follows from the observation that the 
spectral functions as measured in ARPES are always expressed in the plane wave 
basis of the photon, thus in a magnetic field, when the momentum is not a good 
quantum number anymore, it becomes impossible to perform the photoemission 
spectroscopy. 

In order to compute the spectral function, we have to choose a particular fermionic 
plane wave as a probe. Since the separation of variables is valid throughout the bulk, 
the basis transformation can be performed at every constant r-slice. This means that 
only the x and y coordinates have to be taken into account (the plane wave probe 
lives only at the CFT side of the duality). We take a plane wave propagating in the 
+x direction with spin up along the r-axis. In its rest frame such a particle can be 
described by 



Near the boundary (at r^ -^ °°) we can rescale our solutions of the Dirac equation, 
detailes can be found in |fT9l : 



Fi = 



v-^f(.,)c/^'(i)y 






(42) 



with rescaled x defined in |[T9 l. This representation is useful since we calculate the 
components t,±{ri,) related to the retarded Green's function in our numerics (we 
keep the notation of Q). 

Let ffj and ffi be the CFT operators dual to respectively Fi and //, and c|, c^t 
be the creation and annihilation operators for the plane wave state fprobe- Since the 
states F and F form a complete set in the bulk, we can write 



]{(0) 






f/;,c/; 






E 



u;&j{co)+u;^j{co) 



(43) 
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where the overlap coefficients t//((o) are given by the inner product between fprobe 
andF: 

U,ip,) = j dxF} ir'%,^, = -jdxe-P^^^+{rt) (c/''^(i) - C/"^(i)) , (44) 

with F = F^r^, and similar expression for Ui involving ^-(r^). The constants f// 
can be calculated analytically using the numerical value of ^±{ri,), and by noting 
that the Hermite functions are eigenfunctions of the Fourier transform. We are in- 
terested in the retarded Green's function, defined as 

G%{co,p) = -ild''dte'""-'P-^e{t)G%{t,x) 

G%{t,x) ^ {0\[ffiit,x),^,iO,0)]\0) 

where Gff is the retarded Green's function for the operator S'. 

Exploiting the orthogonality of the spinors created by ^ and ^^ and using eq. 
(|43] |. the Green's function in the plane wave basis can be written as 

Gt(ft),P.O = Ltr(^^) (f/*,t7*)G«= (|C/,(/,.)|2G^,(a),0 + |f//(p.)l'G^,(«,0 

(46) 
In practice, we cannot perform the sum in eq. (l46t all the way to infinity, so we have 
to introduce a cutoff Landau level /cut- In most cases we are able to make /cut large 
enough that the behavior of the spectral function is clear 

Using the above formalism, we have produced spectral functions for two dif- 
ferent conformal dimensions and fixed chemical potential and magnetic field (Fig. 
1). Using the plane wave basis allows us to directly detect the Landau levels. The 
unit used for plotting the spectra (here and later on in the paper) is the effective 
temperature Tgs ||5]: 

This unit interpolates between ji at T /jj. = and T and is of or T /ji -^ °o, and is 
convenient for the reason that the relevant quantities (e. g., Fermi momentum) are 
of order unity for any value of /x and h. 



4.2 Magnetic crossover and disappearance of the quasiparticles 

Theoretically, it is more convenient to consider the spectral functions in the Landau 
level basis. For definiteness let us pick a fixed conformal dimension 4 = 1 which 
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Fig. 1 Two examples of spectral functions in the plane wave basis for fi/T = 50 and h/T = 1. The 
conformal dimension is Zi = 5/4 (left) and A = 3/2 (right). Frequency is in the units of effective 
temperature Teff. The plane wave momentum is chosen to be ^' = 1. Despite the convolution of 
many Landau levels, the presence of the discrete levels is obvious. 



corresponds to m = — ^. In the limit of weak magnetic fields, h/T -^ 0, we should 
reproduce the results that were found in |J5]- 

In Fig. (2. A) we indeed see that the spectral function, corresponding to a low 
value of pi/T, behaves as expected for a nearly conformal system. The spectral 
function is approximately symmetric about o = 0, it vanishes for \(o\ < A:, up to a 
small residual tail due to finite temperature, and for \(o\ ^ ^ it scales as oP''" . 

In Fig. (2.B), which corresponds to a high value of ji/T, we see the emergence of 
a sharp quasiparticle peak. This peak becomes the sharpest when the Landau level 
/ corresponding to an effective momentum k^g = ^y2\qh\l coincides with the Fermi 
momentum kp. The peaks also broaden out when k^^ff moves away from kp- A more 
complete view of the Landau quantization in the quasiparticle regime is given in 
Fig. 14.21 where we plot the dispersion relation (co-k map). Both the sharp peaks and 
the Landau levels can be visually identified. 

Collectively, the spectra in Fig. 2 show that conformality is only broken by the 
chemical potential ji and not by the magnetic field. Naively, the magnetic field in- 
troduces a new scale in the system. However, this scale is absent from the spectral 
functions, visually validating the the discussion in the previous section that the scale 
h can be removed by a rescaling of the temperature and chemical potential. 

One thus concludes that there is some value /i^. of the magnetic field, depending 
on n/T, such that the spectral function loses its quasiparticle peaks and displays 
near-conformal behavior for h > h[.. The nature of the transition and the underlying 
mechanism depends on the parameters {pLq,T,A). One mechanism, obvious from 
the rescaling in eq. (|29] l, is the reduction of the effective coupling q as h increases. 
This will make the influence of the scalar potential Aq negligible and push the system 
back toward conformality. Generically, the spectral function shows no sharp change 
but is more indicative of a crossover 

A more interesting phenomenon is the disappearance of coherent quasiparticles 
at high effective chemical potentials. For the special case m = 0, we can go be- 
yond numerics and study this transition analytically, combining the exact T = 
solution found in 1221 and the mapping dSOl l. In the next section, we will show that 
the transition is controlled by the change in the dispersion of the quasiparticle and 
corresponds to a sharp phase transition. Increasing the magnetic field leads to a de- 
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■ (D) - 




Fig. 2 Some typical examples of spectral functions A(tt), ^eff) vs. CO in the Landau basis, ^eff 

5 
4 



i and the bottom four 



y/2\gh\n. The top four correspond to a conformal dimension A 
to 4 = 5 (m = 0). In each plot we show different Landau levels, labelled by index n, as a function 
of ^/T and h/T. The ratios take values {jl/T,h/T) = (1, 1), (50, 1), (1,50), (50,50) from left to 
right. Confomial case can be identified when /i/r is small regardless of h/T (plots in the left 
panel). Nearly confomial behavior is seen when both ^i/T and h/T are large. This confimis our 
analytic result that the behavior of the system is primarily governed by ^. Departure from the 
conformality and sharp quasiparticle peaks are seen when li/T is large and h/T is small in (2.B) 
and (2.F). Multiple quasiparticle peaks arise whenever ^eff = ^'f ■ This suggests the existence of a 
critical magnetic field, beyond which the quasiparticle description becomes invalid and the system 
exhibits a conformal-like behavior As before, the frequency ft) is in units of Teff . 
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-3.0 
-3.5 



(A) 





(B) 



Effective momentum 



Fig. 3 Dispersion relation a vs. ^eff for ^/T = 50, h/T = 1 and 4 = | ('« = ~\)- The spectral 
function A(fi),i'eff) is displayed as a density plot. (A) On a large energy and momentum scale, we 
clearly sees that the peaks disperse almost linearly (ft) r; Vfk), indicating that we are in the stable 
quasiparticle regime. (B) A zoom-in near the location of the Fermi surface shows clear Landau 
quantization. 



crease in phenomenological control parameter v^^ . This can give rise to a transition 
to a non-Fermi liquid when v^^ < 1 /2, and finally to the conformal regime at h = h'^ 
when Vjip ~ and the Fermi surface vanishes. 



4.3 Density of states 



As argued at the beginning of this section, the spectral function can look quite dif- 
ferent depending on the particular basis chosen. Though the spectral function is an 
attractive quantity to consider due to connection with ARPES experiments, we will 
also direct our attention to basis-independent and manifestly gauge invariant quan- 
tities. One of them is the density of states (DOS), defined by 



D(to) = ^A(a),/), 



(48) 



where the usual integral over the momentum is replaced by a sum since only discrete 
values of the momentum are allowed. 

In Fig. 4, we plot the density of states for two systems. We clearly see the Landau 
splitting of the Fermi surface. A peculiar feature of these plots is that the DOS seems 
to grow for negative values of (O. This, however, is an artefact of our calculation. 
Each individual spectrum in the sum eq. ( |48T l has a finite tail that scales as o)^'" for 
large o, so each term has a finite contribution for large values of o. When the full 
sum is performed, this fact implies that lim £)(©) -^ °°. The relevant information 

on the density of states can be obtained by regularizing the sum, which in practice is 
done by summing over a finite number of terms only, and then considering the peaks 
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(A) * 




Fig. 4 Density of states Z)(ffl) for m = - i and (A) n/T = 50,h/T ^ I, and (B) /i/r = 1, h/T = 1 . 
Sharp quasiparticle peaks from the splitting of the Fermi surface are clearly visible in (A). The case 
(B) shows square-root level spacing characteristic of a (nearly) Lorentz invariant spectrum such as 
that of graphene. 



that lie on top of the resuhing finite-sized envelope. The physical point in Fig. 4A 
is the linear spacing of Landau levels, corresponding to a non-relativistic system at 
finite density. This is to be contrasted with Fig. 4B where the level spacing behaves 
as tx: Vh, appropriate for a Lorentz invariant system and realized in graphene | 



5 Fermi level structure at zero temperature 

In this section, we solve the Dirac equation in the magnetic field for the special 
case m = (4 = |). Although there are no additional symmetries in this case, it 
is possible to get an analytic solution. Using this solution, we obtain Fermi level 
parameters such as kp and vp and consider the process of filling the Landau levels 
as the magnetic field is varied. 



5.1 Dirac equation with m = 

In the case m ~ 0, it is convenient to solve the Dirac equation including the spin 
connection (see details in |fT9l ) rather than scaling it out: 



/grr 



-.a^d,-- 



1 



<7 -zfOxfx 



a^{(o + qA,) 






1 



= 0, 



(49) 



where A,/ = y^2|^/i|/ are the energies of the Landau levels Z = 0, 1 , . . ., g,/ = gxx = ^vv 
At{r) is given by eq. Q, and the gamma matrices are defined in l fT9] . In this basis 
the two components i/^i and 1/2 decouple. Therefore, in what follows we solve for 
the first component only (we omit index 1). Substituting the spin connection, we 



Holographic description of strongly coirelated electrons in external magnetic fields 19 

have ES): 



'%/7_n 1 _3,„ , .^ _1'-V7,, , '•/' 



a'dr-^a'ico + qA,)-a'^(3 + ^)-Xijyif = 0, 



R2 
with \j/ = {yi,y2)- It is convenient to change to the basis 

Vi\_M -i\fy^ 



(50) 



~ ; 1 • 1 ; 1 ;' (51) 

which diagonahzes the system into a second order differential equation for each 
component. We introduce the dimensionless variables as in eqs. ( I15H17I ). and make 
a change of the dimensionless radial variable: 

(52) 



1-z' 



with the horizon now being at z = 0, and the conformal boundary at z = 1 . Perform- 
ing these transformations in eq. (l50t . the second order differential equations for _yi 
reads 

+ y{{(o + qfiz)±'-^)^-iq^-xAyi=Q, (53) 

The second component y2 obeys the same equation with /i h^ — /i . 
At r = 0, 



f = 3zHz-zo){z-zo), zo = ^(4 + /\/2). 



(54) 



The solution of this fermion system at zero magnetic field and zero temperature 
T ~ has been found in 1221 . To solve eq. ( f53] l. we use the mapping to a zero 
magnetic field system eq. ( |29] l. The combination jj.^ = jj.q at non-zero h maps to 
;U^ eff = jJ-efSleS ^t zero h as follows: 



/ tt2 / iJ-2 

M./ ^ ?W 1 - q2^h^ ■ 8F\/Q-+H^ = V3qgF\Jl - — = jU,,,eff (55) 

where at T = we used Q^ +H^ ~ 3. We solve eq. ( l53T l for zero modes, i. e. O) = 0, 
and at the Fermi surface X ~k, and implement eq. ( fSSl l. 
Near the horizon (z = 0, / = 6z^), we have 



„2-// , n -/ , / 3 , (M<?.e// 



,2 



6z2y'i';2 + 12zy'i;2 + i^+ ^ '"^^' - 4 bi;2 - 0, (56) 
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which gives the following behavior: 

yi;2-Z"2*^*, (57) 

with the scaling exponent V following from eq. (|32] i: 

v = iy6^M%^, (58) 



at the momentum k. Using Maple, we find the zero mode solution of eq. (I53t with a 
regular behavior z" 2+^ at the horizon 112211181 : 



X 2^1 I - + v-^M,.ff,v + /^:^,l+2v,^-r^ ) , (59) 





and 




X 2^1 ( ^ + V + -^M.;.eff,V~i^^,l+2v " ' "~_ ), (60) 

y2 3 ' 6 3zo(z-zo) ' 

where 2F1 is the hypergeometric function and A^i ,N2 are normalization factors. Since 
normalization factors are constants, we find their relative weight by substituting 
solutions given in eq. i5% back into the first order differential equations at z ~ 0, 

N2 V6k \zoJ 

The same relations are obtained when calculations are done for any z- The second 
solution iilj, with behavior z^ 2^^ at the horizon, is obtained by replacing v -^ —V 
in eq.dlSl). ' 



To get insight into the zero-mode solution ( |59] |. we plot the radial profile for the 
density function i/A^"'^!^'''' for different magnetic fields in Fig. (|5]l. The momentum 
chosen is the Fermi momentum of the first Fermi surface (see the next section). 
The curves are normalized to have the same maxima. Magnetic field is increased 
from right to left. At small magnetic field, the zero modes are supported away from 
the horizon, while at large magnetic field, the zero modes are supported near the 
horizon. This means that at large magnetic field the influence of the black hole to 
the Fermi level structure becomes more important. 
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Fig. 5 Density of the zero mode i/^^ i//" vs. the radial coordinate z (the horizon is at z = and the 
boundary is at j = 1) for different values of the magnetic field /; for the first (with the largest root 



for kf) Fermi surface. We set gf = 1 (/j — > //) and q= -ji. (flci^eff "*■ 15a/ 1 — ^)- From right to left 
the values of the magnetic field are H = {0, 1.40, 1.50, 1.60, 1.63, 1.65, 1.68}. The amplitudes of 
the curves are normalized to unity. At weak magnetic fields, the wave function is supported away 
from the horizon while at strong fields it is supported near the horizon. 



5.2 Magnetic effects on the Fermi momentum and Fermi velocity 



In the presence of a magnetic field there is only a true pole in the Green's function 
whenever the Landau level crosses the Fermi energy lEjj 



2l\qh\ = kp. 



(62) 



As shown in Fig. 2, whenever the equation ( 162b is satisfied the spectral function 
A{(o) has a (sharp) peak. This is not surprising since quasiparticles can be easily 
excited from the Fermi surface. From eq. (l62i , the spectral function A (o) and the 
density of states on the Fermi surface D{(o) are periodic in l with the period 



Ml 



2nq 

a7' 



(63) 



where Ap ~ Ttkj- is the area of the Fermi surface 11231 . This is a manifestation of 
the de Haas-van Alphen quantum oscillations. At T = 0, the electronic properties 
of metals depend on the density of states on the Fermi surface. Therefore, an oscil- 
latory behavior as a function of magnetic field should appear in any quantity that 
depends on the density of states on the Fermi energy. Magnetic susceptibility 1231 
and magnetization together with the superconducting gap ll24ll have been shown 
to exhibit quantum oscillations. Every Landau level contributes an oscillating term 
and the period of the /-th level oscillation is determined by the value of the mag- 
netic field h that satisfies eq. (|62] | for the given value of kp. Quantum oscillations 
(and the quantum Hall effect which we consider later in the paper) are examples of 
phenomena in which Landau level physics reveals the presence of the Fermi sur- 
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face. The superconducting gap found in the quark matter in magnetic fields Il24l is 
another evidence for the existence of the (highly degenerate) Fermi surface and the 
corresponding Fermi momentum. 

Generally, a Fermi surface controls the occupation of energy levels in the sys- 
tem: the energy levels below the Fermi surface are filled and those above are empty 
(or non-existent). Here, however, the association to the Fermi momentum can be 
obscured by the fact that the fermions form highly degenerate Landau levels. Thus, 
in two dimensions, in the presence of the magnetic field the corresponding effective 
Fermi surface is given by a single point in the phase space, that is determined by iif, 
the Landau index of the highest occupied level, i.e., the highest Landau level below 
the chemical potentialQ Increasing the magnetic field. Landau levels 'move up' in 
the phase space leaving only the lower levels occupied, so that the effective Fermi 
momentum scales roughly (excluding interactions) as a square root of the magnetic 
field, kp ~ y/np ^ k"p"' \'i^^TjKiwx- High magnetic fields drive the effective den- 
sity of the charge carriers down, approaching the limit when the Fermi momentum 
coincides with the lowest Landau level. 

Many phenomena observed in the paper can thus be qualitatively explained by 
Landau quantization. As discussed before, the notion of the Fermi momentum is lost 
at very high magnetic fields. In what follows, the quantitative Fermi level structure 
at zero temperature, described by kp and vp values, is obtained as a function of the 
magnetic field using the solution of the Dirac equation given by eqs. ( 159160b . As 
in ifTTI . we neglect first the discrete nature of the Fermi momentum and velocity in 
order to obtain general understanding. Upon taking the quantization into account, 
the smooth curves become combinations of step functions following the same trend 
as the smooth curves (without quantization). While usually the grand canonical en- 
semble is used, where the fixed chemical potential controls the occupation of the 
Landau levels ifZTl , in our setup, the Fermi momentum is allowed to change as the 
magnetic field is varied, while we keep track of the IR conformal dimension v. 

The Fermi momentum is defined by the matching between IR and UV physics 
Q, therefore it is enough to know the solution at o = 0, where the matching is 
performed. To obtain the Fermi momentum, we require that the zero mode solution 
is regular at the horizon (i/A^"' ^ ^^i+^) and normalizable at the boundary. At the 
boundary z ^ 1, the wave function behaves as 

«(l-^)^-'"(j)+Kl-^)^+"'(i)- (64) 

To require it to be normalizable is to set the first term a = 0; the wave function at 
Z ~ 1 is then 

V,«')^(l-z)l+"Yj). (65) 



^ We would like to thank Igor Shovkovy for clarifying the issue with the Fermi momentum in the 
presence of the magnetic field. 
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Fig. 6 Effective momentum kgff vs. the magnetic field h ^ H (we set gp = \, q = ^). As we 
increase magnetic field the Fermi surface shrinks. Smooth solid curves represent situation as if 
momentum is a continuous parameter (for convenience), stepwise solid functions are the real 
Fermi momenta which are discretized due to the Landau level quantization: kf — > yj2\qh\l with 
/ = 1,2, .. . where yj2\qh\l are Landau levels given by dotted lines (only positive discrete kf are 
shown). At a given h there are multiple Fermi surfaces. From right to left are the first, second etc. 
Fermi surfaces. The dashed-dotted line is Vj.j, = where kf is terminated. Positive and negative 
kgff correspond to Fermi surfaces in two components of the Green's function. 




Fig. 7 Landau level numbers n coresponding to the quantized Fermi momenta vs. the magnetic 
field h ^ H for the three Fermi surfaces with positive kf. We sei gf = \, q= M;. From right to 
left are the first, second and third Fermi surfaces. 



Eq. ( l65T l leads to the condition limj.^i (z — 1) ^'^(3'2 + 'Ji ) = 0' which, together 
with eq. (|59] l. gives the following equation for the Fermi momentum as function of 
the magnetic field 12211181 
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2Fl(l + v + %^,i + v-^%£t,l+2v,f(l-/^^)) _ 6v-,X.// 



2Fi(v + %^,i + v-^%^,l+2v,f(l-/V2)) ^f(-2/ + V2)' 



(66) 



with V = Vki, given by eq. (1581 1. Using Mathematica to evaluate the hypergeometric 
functions, we numerically solve the equation for the Fermi surface, which gives 
effective momentum as if it were continuous, i.e. when quantization is neglected. 
The solutions of eq. ( |66] | are given in Fig. (|6]l. There are multiple Fermi surfaces for 
a given magnetic field h. Here and in all other plots we choose gF ~ I, therefore 
h ^)^ H, and ^ = -!4. In Fig.®, positive and negative kp correspond to the Fermi 
surfaces in the Green's functions G\ and G2. The relation between two components 
is G2{(0,k) = Gi (o), —k) [6], therefore Fig.© is not symmetric with respect to the 
X-axis. Effective momenta terminate at the dashed line Vi^p = 0. Taking into account 
Landau quantization of kf — > ^/2\qh\l with Z = 1 , 2 . . ., the plot consists of stepwise 
functions tracing the existing curves (we depict only positive kf). Indeed Landau 
quantiaztion can be also seen from the dispersion relation at Fig. (3), where only 
discrete values of effective momentum are allowed and the Fermi surface has been 
chopped up as a result of it Fig. (3B). 

Our findings agree with the results for the (largest) Fermi momentum in a three- 
dimensional magnetic system considered in 1281 . compare the stepwise dependence 
kpih) with Fig.(5) in Mi- 

In Fig.©, the Landau level index / is obtained from A:^ (/i) ~ y^2|^/j]I where 
kf{h) is a numerical solution of eq. ( |66] |. Only those Landau levels which are below 
the Fermi surface are filled. In Fig.(|6}, as we decrease magnetic field first nothing 
happens until the next Landau level crosses the Fermi surface which corresponds to a 
jump up to the next step. Therefore, at strong magnetic fields, fewer states contribute 
to transport properties and the lowest Landau level becomes more important (see the 
next section). At weak magnetic fields, the sum over many Landau levels has to be 
taken, ending with the continuous limit as h -^ 0, when quantization can be ignored. 

In Fig. ([8]), we show the IR conformal dimension as a function of the magnetic 
field. We have used the numerical solution for kf. Fermi liquid regime takes place 
at magnetic fields h < he, while non-Fermi liquids exist in a narrow band at h^ < 
h < h'e, and at h[ the system becomes near-conformal. 

In this figure we observe the pathway of the possible phase transition exhibited by 
the Fermi surface (ignoring Landau quantization): it can vanish at the line v^^ = 0, 
undergoing a crossover to the conformal regime, or cross the line v^p = 1/2 and go 
through a non-Fermi liquid regime, and subsequently cross to the conformal phase. 
Note that the primary Fermi surface with the highest kp and v^^^ seems to directly 
cross over to conformality, while the other Fermi surfaces first exhibit a "strange 
metal" phase transition. Therefore, all the Fermi momenta with v^^ > contribute 
to the transport coefficients of the theory. In particular, at high magnetic fields when 
for the first (largest) Fermi surface kp is nonzero but small, the lowest Landau 
level « = becomes increasingly important contributing to the transport with half 
degeneracy factor as compared to the higher Landau levels. 
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Fig. 8 Left panel. The IR conformal dimension v = Vj^ calculated at the Fermi momentum vs. 
the magnetic field /i — > // (we set gf =1, q = -^). Calculations are done for the first Fermi surface. 
Dashed line is for v = ^ (at H^ = 1.70), which is the border between the Fermi liquids v > ^ 
and non-Fermi liquids v < ^. Right panel. Phase diagram in terms of the chemical potential and 
the magnetic field /i" + h^ = T> (in dimensionless variables h = gfH, fx = gpQ', we set gf = 1). 
Fermi liquids are above the dashed line {H < H^) and non-Fermi liquids are below the dashed line 
(H > H,). 



In Fig.|9l we plot the Fermi momentum kp as a function of the magnetic field for 
the first Fermi surface (the largest root of eq. (l6&i ). Quantization is neglected here. 
At the left panel, the relatively small region between the dashed lines corresponds to 
non-Fermi liquids < v < j . At large magnetic field, the physics of the Fermi sur- 
face is captured by the near horizon region (see also Fig. (|5)) which is AdS2 x R^. 
At the maximum magnetic field, H,„ia — \/3 w 1 .73, when the black hole becomes 
pure magnetically charged, the Fermi momentum vanishes when it crosses the line 
Vkf = 0. This only happens for the first Fermi surface. For the higher Fermi sur- 
faces the Fermi momenta terminate at the line v^^^ = 0, Fig. (|6]l. Note the Fermi 
momentum for the first Fermi surface can be almost fully described by a function 



H^ 



It is tempting to view the behavior kp ^ \/H,nax — H as a phase 
transition in the system although it strictly follows from the linear scaling for H = 
by using the mapping ( l29l l. (Note that also /i = gfQ = g/r\/3 —H^-) Taking into 
account the discretization of kp, the plot will consist of an array of step functions 
tracing the existing curve. Our findings agree with the results for the Fermi mo- 
mentum in a three dimensional magnetic system considered in [|28l, compare with 
Fig. (5) there. 

The Fermi velocity given in eq. dZTl i is defined by the UV physics; therefore 
solutions at non-zero O) are required. The Fermi velocity is extracted from matching 
two solutions in the inner and outer regions at the horizon. The Fermi velocity as 
function of the magnetic field for v > i is 1221 [TSll 



VF = Y\J "^^^^ 



„^(0)t^(0) 



-1 



lim 



|~(0) 



-iyf? 



(1 
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Fig. 9 Fermi momentum ^^ vs. the magnetic field /i — > // (we set g^ = l,q= -^ ) for the first Fermi 
surface. Left panel. The inner (closer to x-axis) dashed line is Vj.^ = and the outer dashed line is 
Vkp = J, the region between these lines corresponds to non-Fermi liquids < v^^- < j. The dashed- 
dotted line is for the first Landau level k\ = s/lqH. The first Fermi surface hits the border-line 
between a Fermi and non-Fermi liquids v = ^ at Z/^. ks 1.70, and it vanishes at H^ax ~ V3 = 1.73. 
Right panel. Circles are the data points for the Fermi momentum calculated analytically, solid line 

12.96. 



(67) 



is a fit function jf;'"^!- 


^ with ^•™« = 12 


hi 


= lim 





where the zero mode wavefunction is taken at kf eq.i 

We plot the Fermi velocity for several Fermi surfaces in Fig. [TO] Quantization 
is neglected here. The Fermi velocity is shown for v > 5. It is interesting that the 
Fermi velocity vanishes when the IR conformal dimension is v^.^ = j. Formally, 
it follows from the fact that v/r ^ (2v — 1) Q. The first Fermi surface is at the 
far right. Positive and negative vp correspond to the Fermi surfaces in the Green's 
functions Gi and G2, respectively. The Fermi velocity Vf has the same sign as the 
Fermi momentum kp. At small magnetic field values, the Fermi velocity is very 
weakly dependent on H and it is close to the speed of light; at large magnetic field 
values, the Fermi velocity rapidly decreases and vanishes (at He = 1.70 for the first 
Fermi surface). Geometrically, this means that with increasing magnetic field the 
zero mode wavefunction is supported near the black hole horizon fig.®, where the 
gravitational redshift reduces the local speed of light as compared to the boundary 
value. It was also observed in Il22l l7l at small fermion charge values. 



6 Hall and longitudinal conductivities 



In this section, we calculate the contributions to Hall (7^. and the longitudinal (7vx 
conductivities directly in the boundary theory. This should be contrasted with the 
standard holographic approach, where calculations are performed in the (bulk) grav- 
ity theory and then translated to the boundary field theory using the AdS/CFT dictio- 
nary. Specifically, the conductivity tensor has been obtained in y_0'| by calculating 
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Fig. 10 Fenni velocity Vf vs. the magnetic field h ^ H (we set gf = 1, ly = -^) for the regime 

of Fermi liquids V > j- Fermi velocity vanishes at Vkf = 5 (x-axis). For the first Fermi surface, 
the top curve, fermi velocity vanishes at H^ Ki 1.70. The region H < He corresponds to the Fermi 
liquids and quasiparticle description. The multiple lines are for various Fermi surfaces in ascending 
order, with the first Fermi surface on the light. The Fermi velocity Vf has the same sign as the 
Fermi momentum kf. As above, positive and negative v^ correspond to Fermi surfaces in the two 
components of the Green's function. 



the on-shell renormalized action for the gauge field on the gravity side and using 
the gauge/gravity duality Am -^ j/i to extract the R charge current-current correlator 
at the boundary. Here, the Kubo formula involving the current-current correlator is 
used directly by utilizing the fermion Green's functions extracted from holography 
in Q. Therefore, the conductivity is obtained for the charge carriers described by 
the fermionic operators of the boundary field theory. 

The use of the conventional Kubo formulo to extract the contribution to the trans- 
port due to fermions is validated in that it also follows from a direct AdS/CFT com- 
putation of the one-loop correction to the on-shell renormalized AdS action ifTSll . 
We study in particular stable quasiparticles with V > 5 and at zero temperature. 
This regime effectively reduces to the clean limit where the imaginary part of the 
self-energy vanishes ImE -^ 0. We use the gravity-"dressed" fermion propagator 
from eq. jTU and to make the calculations complete, the "dressed" vertex is neces- 
sary, to satisfy the Ward identities. As was argued in IfTSl , the boundary vertex which 
is obtained from the bulk calculations can be approximated by a constant in the low 
temperature limit. Also, according to ||26) , the vertex only contains singularities of 
the product of the Green's functions. Therefore, dressing the vertex will not change 
the dependence of the DC conductivity on the magnetic field 1261 . In addition, the 
zero magnetic field limit of the formulae for conductivity obtained from holography 
IfTSl and from direct boundary calculations IJTSl are identical. 
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6.1 Integer quantum Hall effect 

Let us start from the "dressed" retarded and advanced fermion propagators Q: Gr 
is given by eq. (|27] | and Ga = G^. To perform the Matsubara summation we use the 
spectral representation 

GM.k)=/^^, (68) 

J 2k (0 — l(On 

with the spectral function defined as A{co,k) = —j^lmGR{co,k) — ^■{GR{co,k) — 
GA{(0,k)). Generalizing to a non-zero magnetic field and spinorcase ll27l . the spec- 
tral function ll29l is 



1 _*' 
A{co,k) = -e W (69) 



where Ef = vf^f is the Fermi energy, £/ ~ VF\/2\qli\l is the energy of the Landau 

level, /(k) = P^Li^Yji) — P+Li_\{tjt) with spin projection operators P± = (1 ± 

iy^'f')/!, we take c = 1, the generalized Laguerre polynomials are L"(z) and by 
definition L„ (z) = L[J (z), (we omit the vector part ky, it does not contribute to the DC 
conductivity), all 7's are the standard Dirac matrices, hy, Vf and kf are real constants 
(we keep the same notations for the constants as in Q)- The self-energy L ^ CO ^''f 
contains the real and imaginary parts, E = Ei+ iE2. The imaginary part comes from 
scattering processes of a fermion in the bulk, e. g. from pair creation, and from the 
scattering into the black hole. It is exactly due to inelastic/dissipative processes that 
we are able to obtain finite values for the transport coefficients, otherwise they are 
formally infinite. 

Using the Kubo formula, the DC electrical conductivity tensor is 

Imnff 
(j,/(f2)= hm ^, (70) 

where IIij{iQ,n -^ £2 +/0+) is the retarded current-current correlation function; 
schematically the current density operator is /(t,x) ~ qvpY^a '/fj(T,x)7'\//(j(T,x). 
Neglecting the vertex correction, it is given by 

nij[iQ,n) = q^\^FT £ / T^tr {YG{i(0„,k)y^G{i(0„ + /A„,k)) . (71) 

The sum over the Matsubara frequency is 

tY — 1 _ «(coi)-«(a)2) ^^2) 

„ iO)n — COl iCOn + Him — 0>2 /i2,„ + COi — (02 
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Taking iQ,„ -^ Q + iO^, the polarization operator is now 

d(0id0>2npD{(0i)-npD{0}2) f (fk , ■ ■ . 

^'■^■^^) - -2^^ X2 + a),-a>, I (2^tr(/A(co„k)7M(a,„k)) ,(73) 

where the spectral function A(«,k) is given by eq. ( |69] l and «FD(fi^) is the Fermi- 
Dirac distribution function. Evaluating the traces, we have 

= _JL^^±I1J1A (74) 

TTii 

X Re £ {-!)'+'+' {8ij{5,,k-i+d,^i,k)+ieijSgn{qh){8,,,^i-5i_i,,)} 

l,k=0 

(tanh tanh — ) 

(a„-£,)2 + £|(a),)^ ' ""V I (0)2- £.)2 + 2:1(0)2) + ^' '^ 

with CO2 = coi + Q . We have also introduced (bia = (Oi;2 + £f + ^i(fi)i;2) with e,y 
being the antisymmetric tensor (£12 = 1), and Zi;2(fi)) = i;i;2(ft),^f ). In the mo- 
mentum integral, we use the orthogonality condition for the Laguerre polynomials 
Jq dxt''Li{x)Lk{x) = 5ik. 

From eq. ( |74l i. the term symmetric/antisymmetric with respect to exchange 
(0\ -^ (i>2 contributes to the diagonal/off-dialgonal component of the conductivity 
(note the antisymmetric term npoCci) — nFD(fi)2))- The longitudinal and Hall DC 
conductivities {Q. — > 0) are thus 



2q\hxVFY\qh\ n dco r|(«) 
Irf 7-00 27: cosh^^ 



a„ = -rLSl^^LUim I l^J=2^ (75) 



/= 



'o\{d)-Eiy + E^i(o) 'J \{cb-Ei+iy+E^{co) 



xI( .^,_J^.2,.., +(^/^-^/))( ,^,_.,..:2^.2,.., +(^m^-£m)), 



q^{hiVF)^sgn{qh) 

(y„ ^ V/,, 

n 

vi,=2[ ^tanh:^i:2(co)f a,f--— -^-— 27— + (£/^-£/) 
J-ooIn 2T f^Q \{co-Eiy + E^[co) 



where cb = CO + £f + Ei {(o)). The filling factor V/, is proportional to the density 
of carriers: |v/,| = T^m;-^n (see derivation in lfT9l ). The degeneracy factor of the 
Landau levels is a/: Oq = 1 for the lowest Landau level and «/ = 2 for / = 1,2. . .. 
Substituting the filling factor V/, back to eq. ( |76] l, the Hall conductivity can be written 

^' P 

(Tn. - 7, (77) 

h 



(76) 
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where p is the charge density in the boundary theory, and both the charge q and 
the magnetic field h carry a sign (the prefactor {—h\Vf) comes from the normal- 
ization choice in the fermion propagator eqs. (I27l69t as given in [|7|, which can be 
regarded as a factor contributing to the effective charge and is not important for 
further considerations). The Hall conductivity eq. (l77t has been obtained using the 
AdS/CFT duality for the Lorentz invariant 2 + 1 -dimensional boundary field theo- 
ries in ifTol . We recover this formula because in our case the translational invariance 
is maintained in the x and y directions of the boundary theory. 

Low frequencies give the main contribution in the integrand of eq. (|76] |. Since 
the self-energy satisfies Z\ {(o) ^ ^^2(0) ^ C^*' and we consider the regime V > ^, 
we have Zi ~ Z2 ^" at © ~ (self-energy goes to zero faster than the (O term). 
Therefore, only the simple poles in the upper half-plane oq = — Ef ± £/ + Zj + /Z2 
contribute to the conductivity where Zi ~^ Z2 ^ (— Ef ±£/)^*' are small. The same 
logic of calculation has been used in ll27l . We obtain for the longitudinal and Hall 
conductivities 

2,^(/..v.)^Z2 / 1 £^^l+coshfcoshf^Y 

nT \^l+coshf ,ti (coshf +cosh^)2y' 

a^, ^ \^ ^ 2 I tanh— + £(tanh^^ +tanh^^) 1 , (79) 

where the Fermi energy is £f ~ Vfkp and the energy of the Landau level is £/ = 
VF\/2\qh\l. Similar expressions were obtained in lIZTl . However, in our case the 
filling of the Landau levels is controlled by the magnetic field h through the field- 
dependent Fermi energy VF{h)kF{h) instead of the chemical potential /i. 

At r = 0, cosh^ — > ^eT" and tanh^ = 1 — 2n-pT){(o) -^ sgnco. Therefore the 
longitudinal and Hall conductivities are 

2q^{h,VFfZ2^,^ 2q^{h,VpfZ2 <, 



(y_„ = 



q^{hiVffsgn{qh) 
n 

q'^{hiVffsgn{qh) 

n 



2ll+2^0(ef-£,)J 
x2(l+2n)0(ef-£„)0(£„+i-ef), (81) 



where the Landau level index runs n = 0, 1 , It can be estimated as « = t, , 

when \'F 7^ ([ ] denotes the integer part), with the average spacing between the 
Landau levels given by the Landau energy VF\j2\qh\. Note that £f = Spih)- We 
can see that eq. dsTl i expresses the integer quantum Hall effect (IQHE). At zero 
temperature, as we dial the magnetic field, the Hall conductivity jumps from one 
quantized level to another, forming plateaus given by the filling factor 
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V/, = ±2(l+2n) = ±4(n + i), (82) 

with n = 0, 1, (Compare to the conventional Hall quantization V/, = ±4n, that 

appears in thick graphene). Plateaus of the Hall conductivity atT = follow from 
the stepwise behavior of the charge density p in eq.dTTli: 

p ^ 4(« + i)0(eF -£„)0(£„+i - Sf), (83) 

where n Landau levels are filled and contribute to p. The longitudinal conductivity 
vanishes except precisely at the transition point between the plateaus. In Fig.[TTl we 
plot the longitudinal and Hall conductivities at T = 0, using only the terms after x 
sign in eq. ( |79] ). In the Hall conductivity, plateau transition occurs when the Fermi 
level (in Fig. fTTT i of the first Fermi surface Ep — Vf{h)kf{li) (Fig. |9]l crosses the 
Landau level energy as we vary the magnetic field. By decreasing the magnetic 
field, the plateaus become shorter and increasingly more Landau levels contribute 
to the Hall conductivity. This happens because of two factors: the Fermi level moves 
up and the spacing between the Landau levels becomes smaller This picture does 
not depend on the Fermi velocity as long as it is nonzero. 



6.2 Fractional quantum Hall effect 

In |fT9ll . using the holographic description of fermions, we obtained the filling factor 
at strong magnetic fields 

V/,=±2;, (84) 

where j is the effective Landau level index. Eq.(|84]i expresses the fractional quan- 
tum Hall effect (FQHE). In the quasiparticle picture, the effective index is integer 
7 = 0, 1,2, • • •, but generally it may be fractional. In particular, the filling factors 
V = 2/m where m — 1,2,3, •• • have been proposed by Halperin ll33l for the case 
of bound electron pairs, i.e. 2e-charge bosons. Indeed, QED becomes effectively 
confining in ultraquantum limit at strong magnetic field, and the electron pairing is 
driven by the Landau level quantization and gives rise to 2e bosons. In our holo- 
graphic description, quasiparticles are valid degrees of freedom only for v > 1/2, 
i.e. for weak magnetic field. At strong magnetic field, poles of the fermion propaga- 
tor should be taken into account in calculation of conductivity. This will probably 
result in a fractional filling factor Our pattern for FQHE Fig.lfTTTi resembles the one 
obtained by Kopelevich in Fig. 3 |32il which has been explained using the fractional 
filUng factor of Halperin ll33l . 

The somewhat regular pattern behind the irregular behavior can be understood 
as a consequence of the appearance of a new energy scale: the average distance 
between the Fermi levels. For the case of Fig. [TT] we estimate it to be < Ep — 
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^(m+\) ^^4 9 ^jjjj OT :^ 1 2. The authors of W\ explain the FQHE through the 
opening of a gap in the quasiparticle spectrum, which acts as an order parameter 
related to the particle-hole pairing and is enhanced by the magnetic field (magnetic 
catalysis). Here, the energy gap arises due to the participation of multiple Fermi 
surfaces. 

A pattern for the Hall conductivity that is strikingly similar to Fig.dTTTi arises in 
the AA and AB-stacked bilayer graphene, which has different transport properties 
from the monolayer graphene llSTI . compare with Figs. (2,5) there. It is remarkable 
that the bilayer graphene also exhibits the insulating behavior in a certain parameter 
regime. This agrees with our findings of metal-insulating transition in our system. 
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Fig. 11 Hall conductivity a^. and longitudinal conductivity (T„ vs. the magnetic field h - 
r = (we set gf = I, q = ^). Left panel is for IQHE. Right panel is for FQHE. At strong 
magnetic fields, the Hall conductivity plateau V/, = 4 appears together with plateaus V;, = 2 and 
V;, = 6 in FQHE (detailes are in 1191 ). Irregular pattern in the length of the plateaus for FQHE is 
observed in experiments on thin films of graphite at strong magnetic fields |25| . 



7 Conclusions 



We have studied strongly coupled electron systems in the magnetic field focussing 
on the Fermi level structure, using the AdS/CFT correspondence. These systems are 
dual to Dirac fermions placed in the background of the electrically and magnetically 
charged AdS-Reissner-Nordstrom black hole. At strong magnetic fields the dual 
system "lives" near the black hole horizon, which substantially modifies the Fermi 
level structure. As we dial the magnetic field higher, the system exhibits the non- 
Fermi liquid behavior and then crosses back to the conformal regime. In our analysis 
we have concentrated on the the Fermi liquid regime and obtained the dependence of 
the Fermi momentum kp and Fermi velocity Vf on the magnetic field. Remarkably, 
kp exhibits the square root behavior, with vp staying close to the speed of light in a 
wide range of magnetic fields, while it rapidly vanishes at a critical magnetic field 
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which is relatively high. Such behavior indicates that the system may have a phase 
transition. 

The magnetic system can be rescaled to a zero-field configuration which is ther- 
modynamically equivalent to the original one. This simple result can actually be 
seen already at the level of field theory: the additional scale brought about by the 
magnetic field does not show up in thermodynamic quantities meaning, in particu- 
lar, that the behavior in the vicinity of quantum critical points is expected to remain 
largely uninfluenced by the magnetic field, retaining its conformal invariance. In the 
light of current condensed matter knowledge, this is surprising and might in fact be 
a good opportunity to test the applicability of the probe limit in the real world: if 
this behavior is not seen, this suggests that one has to include the backxeaction to 
metric to arrive at a realistic description. 

In the field theory frame, we have calculated the DC conductivity using kf and 
vf values extracted from holography. The holographic calculation of conductivity 
that takes into account the fermions corresponds to the corrections of subleading 
order in 1 /N in the field theory and is very involved ifTSi . As we are not interested 
in the vertex renormalization due to gravity (it does not change the magnetic field 
dependence of the conductivity), we have performed our calculations directly in the 
field theory with AdS gravity-dressed fermion propagators. Instead of controlling 
the occupancy of the Landau levels by changing the chemical potential (as is usual 
in non-holographic setups), we have controlled the filling of the Landau levels by 
varying the Fermi energy level through the magnetic field. At zero temperature, we 
have reproduced the integer QHE of the Hall conductivity, which is observed in 
graphene at moderate magnetic fields. While the findings on equilibrium physics 
(Landau quantization, magnetic phase transitions and crossovers) are within expec- 
tations and indeed corroborate the meaningfulness of the AdS/CFT approach as 
compared to the well-known facts, the detection of the QHE is somewhat surpris- 
ing as the spatial boundary effects are ignored in our setup. We plan to address this 
question in further work. 

Interestingly, at large magnetic fields we obtain the correct formula for the fill- 
ing factor characteristic for FQHE. Moreover our pattern for FQHE resembles the 
one obtained in ll32l which has been explained using the fractional filling factor of 
Halperin ll33l . In the quasiparticle picture, which we have used to calculate Hall 
conductivity, the filling factor is integer In our holographic description, quasiparti- 
cles are valid degrees of freedom only at weak magnetic field. At strong magnetic 
field, the system exhibits non-Fermi liquid behaviour In this case, the poles of the 
fermion propagator should be taken into account to calculate the Hall conductivity. 
This can probably result in a fractional filling factor We leave it for future work. 

Notably, the AdS-Reissner-Nordstrom black hole background gives a vanishing 
Fermi velocity at high magnetic fields. It happens at the point when the IR conformal 
dimension of the corresponding field theory is v = ^, which is the borderline be- 
tween the Fermi and non-Fermi liquids. Vanishing Fermi velocity was also observed 
at high enough fermion charge ll22l . As in ||221 . it is explained by the red shift on 
the gravity side, because at strong magnetic fields the fermion wavefunction is sup- 
ported near the black hole horizon modifying substantially the Fermi velocity. In 
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our model, vanishing Fermi velocity leads to zero occupancy of the Landau levels 
by stable quasiparticles that results in vanishing regular Fermi liquid contribution 
to the Hall conductivity and the longitudinal conductivity. The dominant contribu- 
tion to both now comes from the non-Fermi liquid and conformal contributions. 
We associate such change in the behavior of conductivities with a metal-"strange 
metal" phase transition. Experiments on highly oriented pyrolitic graphite support 
the existence of a finite "offset" magnetic field he at T = where the resistivity 
qualitatively changes its behavior ll30l . At T 7^ 0, it has been associated with the 
metal-semiconducting phase transition |f30l. It is worthwhile to study the temper- 
ature dependence of the conductivity in order to understand this phase transition 
better 
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